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Abstract 

We explore the sensitivity of flavour changing b — > s transitions to a (sub-)GeV hidden 
sector with generic couplings to the Standard Model through the Higgs, vector and axion 
portals. The underlying two-body decays of B mesons, B — > X S S and B° — > SS, where S 
denotes a generic new GeV-scale particle, may significantly enhance the yield of monochro- 
matic lepton pairs in the final state via prompt S —> 11 decays. Existing measurements of 
the charged lepton spectrum in neutral-current semileptonic B decays provide bounds on the 
parameters of the light sector that are significantly more stringent than the requirements of 
naturalness. New search modes, such as B — > X s + n(ll) and B° — > n(ll) with n > 2, can 
provide additional sensitivity to scenarios in which both the Higgs and vector portals are 
active, and are accessible to (super-)i? factories and hadron colliders. 



1. Introduction 



The study of B mesons at the 5-factories, BaBar pQ and Belle [2], and the Tevatron ex- 
periments [3111] has significantly advanced the precision with which various Standard Model 
(SM) parameters are known, and consequently has placed stringent constraints on models of 
new physics affecting quark flavour [5]. The prevailing view is that such new physics must 
reside at or above the electroweak scale, manifesting at low energies in modifications to the 
Wilson coefficients of effective flavour-changing operators that arise once the heavy degrees 
of freedom are integrated out. Experimental precision, and the ability to make accurate 
SM predictions, are thus the controlling factors in probing weak-scale new physics through 
precision flavour observables. 

While new states charged under the SM are generically required to be rather heavy, light 
(sub-)GeV mass states in a hidden sector, neutral under the SM gauge group, can peacefully 
co-exist with the SM, evading precision flavour and electroweak constraints. Such hidden 
sectors may be weakly coupled to the SM in various ways, and are often best probed via 
experiments at the luminosity frontier. In particular, precision studies of rare SM decays 
can provide impressive sensitivity to these sectors, opening the possibility for novel decay 
channels not encountered in the SM itself. Indeed, over the years there have been numerous 
searches for rare decays of flavoured mesons to new light states (see e.g. [6] for a subset 
of theoretical ideas). As one notable motivation, these hidden sector states can have a 
significant impact on Higgs decay channels, allowing for a SM-like Higgs with mass well 
below the conventional LEP bound [7]. 

In this paper, we revisit the sensitivity of rare flavour-changing decays from the generic 
standpoint of 'portal' operators [SJ[5], which constitute a systematic way to parametrize 
the allowed couplings of generic neutral states S in a hidden sector to the SM in order of 
increasing canonical operator dimension. In particular, we will be interested in the following 
set of lowest- dimension portals: 

H ] H(AS + XS 2 ) Higgs portal (dim = 3, 4), 

kF Iu F Iu Vector portal (dim = 4), (1) 

YnLHN Neutrino portal (dim = 4), 

fa 1 ' l Pl^l5^pd^a Axion portal (dim = 5). 

Here H is the SM Higgs doublet, F^ u is the hypercharge field strength, L is the left-handed 
lepton doublet, and ip is a generic SM fermion, while S = S, N, A'^ and a denote the fields 
associated with new light states. The purpose of this study is to analyze the feasibility of 
searching for light states coupled to the SM via these portals in B meson decaysQ Specifically, 
we will concentrate on the manifestations of Higgs, vector, and axion portals in b — > s 
transitions with the direct production of one or more exotic states. To be as conservative 
as possible, we shall not assume any direct flavour- violating operators, which in fact is 
automatic for the Higgs and vector portals, but requires an extra assumption for the axion 



1 Renewed interest in the possibility of light hidden sector states coupled to the SM has emerged from 
attempts to link certain unexpected features in the multi-GeV scale cosmic electron and positron spectra to 
the annihilation of galactic dark matter into such light states [iOlITT] . 
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portal. Using the resulting flavour-blind portal operators, we calculate the strength of the 
flavour- changing transitions induced by SM loops. 

A primary feature that we will exploit is that the scalar and axion (i.e. axial-vector) 
portals behave very differently to the conserved vector current portal once dressed by W — 
(u, c, t) loop corrections. Schematically, this difference can be illustrated as follows: 

*7m* — y ( G Fq 2 ) x Hj^sl] tj^t — > {G F m 2 t ) x b L ^^s L . (2) 

While conservation of the vector current (such as the electric charge or baryon number) 
requires the dependence on q 2 < ml, the axial current is not conserved and the vertex 
correction is 0(m 2 /q 2 )-enhanced relative to the vector case. Within the SM, scalar or axial- 
vector currents are associated purely with couplings to the Z boson and the SM Higgs, which 
cannot be produced in on-shell B decays. Thus, having light states in the spectrum with 
(pseudo) scalar or axial- vector couplings can enhance the loop-induced two-body decays of 
the b quark by many orders of magnitude. The enhancement of the loop-induced SM Higgs 
coupling has been known for some time [T2"|lf3]. More recently it has been exploited in 
the context of B meson decays to a pair of light dark matter particles through the Higgs 
portal [H], decays to a singlet scalar mixed with the Higgs [T5], and decays to a light 
pseudoscalar in the NMSSM [IB]. Rare Kaon decays to metastable mediators were considered 
in [T71CE8]. 

We will analyze a number of semi-leptonic and fully leptonic B decay modes opened 
up by portal couplings, which can serve as a powerful probe of new light states. As often 
happens in models of this type with intermediate cascade decays, the increased multiplicity 
of final state leptons implies minimal additional suppression [T8H20] . thus enhancing signal 
over background. Specifically, we calculate B -)■ K(K*)S -> K(K*)ll and B° ->■ SS -> 2(11) 
in the minimal extension of the SM by one real scalar S, and B — > K(K*)a — > K(K*)lT 
in the axion portal model. We will show that the constraints imposed by 5-physics in the 
kinematically accessible range where the leptonic decays of S and a occur within the detector 
are easily the most stringent experimental limits. We also extend our analysis to include 
the vector portal, and in particular the natural combination of Higgs and vector portals, 
and calculate the branching of B — > VK(K*), B° — > VV and B —> h'h'. The final state of 
two Higgs h' bosons of the extra U(l) group may be dominated by eight leptons. The most 
important point of our analysis is to show that multilepton signatures of B meson decays, 
like B° — > fi~ fi + fi~ fi + , can be explored using existing datasets collected at the B factories 
and the Tevatron, providing significant new probes of these models with exotic light neutral 
states. 

The remainder of this paper is structured as follows. In section 2, we analyze rare B 
decays in the minimal extension of the SM by a singlet scalar interacting through the Higgs 
portal, as well as an extension with a pseudoscalar singlet coupled via the axion-portal. 
Section 3 considers rare B decay modes proceeding via a combination of Higgs and vector 
portals, and we present our conclusions in section 4. 
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2. Rare 5-decays through the Higgs and axion portals 

The extension of the SM by a singlet scalar has been considered on numerous occasions, e.g. 
for cosmological applications as a minimal model of dark matter, with stability imposed by 
symmetry [2U[22], or its impact on electroweak baryogenesis or inflation [23]. Novel exper- 
imental signatures, including extra decay channels for the SM Higgs boson, were addressed 
in 0[ig[22j|2lH28]. 

A generic renormalizable scalar potential that includes S self-interactions and couplings 
to the SM via the Higgs portal is given by, 

V = A4S 4 + A 3 S 3 + m 2 S 2 + {AS + XS 2 )(H^H). (3) 

Since we are interested only in the low-energy limit of the theory relevant for B decays, we 
will assume stability of the potential in the S-direction and integrate out the Higgs boson to 
obtain an effective Lagrangian for S (enforcing (S) = by an appropriate shift of the field), 

£.= i(V)"-im|S"-(f + ^)^-^ + -. (4) 

The quantity C m comprises the SM mass terms from electroweak symmetry breaking (i.e. 
C m = mill + •••), and the physical mass ms, mixing angle 9, and self-interaction parameter 
A' are related to the parameters in ([3]). The precise nature of these relations [9 ~ Avjm\ 
etc.) will not be critical to our analysis. However, the technical naturalness of the model 
(jlj) is a valuable criterion to use in setting the characteristic values of 9 and A'. In order to 
shelter a relatively light scalar from large mass corrections induced by electroweak symmetry 
breaking, we take 

< ~ O(10") x , 
m h VI GeV/ 

A' < (16vr 2 m|) 1/2 ~ 0(10 GeV) x (J^) . (5) 

The latter relation follows from the SS loop correction to the mass of the scalar. A larger 
angle 9 and self-interaction parameter A' would require additional tuned cancellations be- 
tween different contributions to ms. The possibility of a stronger coupling to the Higgs 
portal, while keeping S light and avoiding the naturalness constraints, arises in the large 
tan/3 two-Higgs doublet extension of the SM [H] and thus also in the MSSM. 

We will also explore the axion portal, which avoids corrections to the (sub-) GeV mass of 
the pseudoscalar via the dimension-five axial-vector couplings of the form, 

C a = £ ^l.lH- (6) 
SM-V> ^ 

Furthermore, for simplicity, we will neglect the effects of the self-interaction of a, as well as 
couplings to gauge bosons, and assume universal couplings of the pseudoscalar to leptons 
fi and quarks, f q . This automatically protects (EJ) from tree-level flavour changing neutral 
currents (FCNCs). While a UV completion is required for ([6]), we note that in two-Higgs 
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doublet extensions of the SM there also exists the possibility of a renormalizable pseudoscalar 
portal, e.g. iaHiH 2 , which leads to the mixing of a with the pseudoscalar Higgs boson A. 

For both the Higgs and axion portals, on integrating out the H^-top loop, we obtain the 
well-known effective b — s — h and b — s — a vertices, 



3V2G F m 2 t V t *V tb , [OS ,2 a, /l2 , ,A , 

bs = 16tt2 X mb " L r X [~V ~ l 3j ( A uv/ m ?) J + ( h - c 



(7) 



For the scalar S, the Wilson coefficient in ([7]) is one- loop exact in the limit m 2 /M^ — > 0, 
while for the pseudoscatal we retain only the leading log-divergent term proportional to 
mj/rriyy and for consistency assume at least a small hierarchy between the weak scale and 
the UV cutoff, hiAuv/mt ~ 1. We have integrated by parts and used the equations of 
motion for the quark fields in the limit m s = to remove the derivative from the axion field 
in the interaction (J7|). 

The Lagrangian (J7|) immediately leads to the inclusive b quark decay width to S and a, 
but we are more interested in K and K* final states. The QCD matrix elements involved 
in -Brf( u ) to K(K*) transitions have been calculated using light-cone QCD sum rules [30|[3Tj. 
and after a fairly standard calculation, we obtain the following results as functions of mj 



and m a : 



Btb^ks ^ 4 x 1(T 7 x (JL^J JFl{ms)\f s 

B rjw5 ~ 5 x 10- 7 x T 2 K *{m s )\%l s 

6 /lOOTeV fA VY \\ 2 _ 2 . 1/2 

BlB^Ka ^ 5 X 10 X ^ In J ^K{ m a)^Ka 
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100 TeV, Muv\ \ ^2 , ,,3/2 



Br B -tf.« ~ 6xl0- b x^— ln V^fJJ F K^a)\%t a - 

The dependence on the unknown mass parameters resides in the phase space factors, Ay = 
(1 — m^ 2 (mj + rrij) 2 )(l — m^ 2 (mj — mj) 2 ), and the form factors which we have normalized 
to their values at zero momentum transfer |31j . 

_ i 

F K {m) 



2\ ' 



1 - m 2 /(38 GeV 

3.65 2.65 

K * [m) ~ 1 - m 2 /(28 GeV 2 ) ~ 1 - m 2 /(37 GeV 2 ) ' U 

The values of the form factors at q 2 = used in our calculations are /o(0) = 0.33 and 
A (0) = 0.37 [31]. The uncertainty in the form factors is the main source of error for (jSJ), 
argued to be at the O(10 — 15%) level 



2 We thank the authors of Ref. [29] for pointing out the presence of a logarithmic UV divergence in this 
calculation. 



4 



The results in (JS}, combined with the subsequent decay of a or S to dilepton pairs 
close to the interaction point presents an intriguing signal: a monoenergetic lepton pair in 
association with K or K*. The branching ratios Bre^^/j — 4.2+Qg x 10~ 7 and Btb^-k*^ — 
1.031 

a23 x 10 6 have been measured [32H3I] with several hundred decays containing lepton 
pairs distributed over the entire available q 2 range, while a monoenergetic lepton pair can 
be efficiently probed at I?-factories with O(l0~ 8 ) sensitivity [35]. The hadronic decays of 
S and a as well as missing energy signatures from decays outside the detector can also be 
probed, albeit with lesser sensitivity. 

With S and a in the intermediate state, the decay widths and branching ratios to leptons 
are sensitive functions of mass. We follow the standard prescriptions for calculating the total 
widths of a and 5" [36H38] . and the results can be summarized as follows. When only decays 
to leptons are kinematically allowed, the leptonic branching is necessarily close to unityjf] 
while the partial decay width to a lepton pair is given by, 

(Pmfmsf Am 2 \ 3/2 mfm a ( Am 2 \ 1/2 

871V 2 V m S J 87T fl \ m a J 

and is very sensitive to whether the dimuon channel is open. For example, for a 250 MeV 
mass scalar with mixing angle 10~ 3 the lifetime is cr = 2.7 cm, and considering a Lorentz 
boost of 7 ~ tub/ (2ms) ~ 10, this would correspond to a significantly displaced vertex. 

For higher mass scalars, the decay length shrinks while the leptonic branching gets sup- 
pressed, especially near the /o + resonance [38]. In the region near the resonance, we 
base our estimate of the branching on a coupled-channel analysis in the framework of chiral 
perturbation theory, while above the resonance, we use perturbative QCD [36|l38]: 



TV? 

Brs-^/j ~ on* . Ir , , n tTTT , i r, /„ TT7T- for m s <1.5GeV, 



rn 



2/33 + \F n /2m s \ 2 ^ + \F K /2m s \ 2 f3 K 



m l 



where fa — (1 — 4m^/m|) 1/ ' 2 , are the form-factors defined in Ref. [SB], Nf is the number 
of heavy quarks, i.e. three below the charm threshold, and the ellipsis in the second line 
stands for charm and r contributions once the corresponding thresholds are open. We note 
that there is at least a 100% uncertainty in this formula above 1 GeV [36] . 

For the pseudoscalar case, the hadronic width is suppressed by three pion phase space. 
In order to estimate the scaling of the branching ratio with fi and f q we assume that the 
decay to hadrons occurs via mixing with the t] and i]' resonances. Taking a representative 
value of m a = 800 MeV, the mixing with rj' is given by 9 ari ' ~ (f v ' / fq) x y/^ m a/( m2 ,> ~ m a), 
and the hadronic width is approximately I^ad ~ @a V '^v'- Using these results, we obtain the 
following scaling of the leptonic branching fraction: 

Bra ^ A ~ l + 0.3(/,// 9 )2- (12) 



3 Within the axion portal scenario, for certain parameter choices the decay to 77 can be significant and 
may also be a good search mode. 
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0.2 0.5 1.0 2.0 5.0 

m s (GeV) 

Figure 1: Sensitivity of the BaBar/Belle dataset to combined B — > KS and B — > K*S decays in the 
dimuon channel. The region below the dashed line is technically natural as discussed in Eq. (J5J) - 

It is apparent that the resonant enhancement of the hadronic width can significantly exceed 
the naive three-pion continuum result. 

With these estimates in hand, we can predict the observable signal at (super-) B factories. 
Having a typical detector design in mind, we require S or a to decay within a transverse 
distance l m i Q = 25 cm of the beam pipe, and assume ~ 90% angular acceptance. In practise 
this amounts to calculating the following angular integral multiplying Eqs. (jHJ): 

f*- *** sin 9d9 

Br S ( )_^ / ^ 

J 8 mill 

In the limit of a short decay length, the integral is trivially Brs( )_> M £ cos 6* min , and in the 
opposite limit of a very long decay length it is (r 's(a)-*np l mm) x (tt/2 - 9 min )j~^ a y 

Given that the combined BaBar/Belle dataset provides sensitivity to the Kfifl and K*n/2 
branching with a mono-energetic muon pair at the level of O(10~ s ), the significant parameter 
space reach that ensues for the two models is shown in Figs. Q] and El For the scalar singlet 
Higgs portal, Fig. [1] illustrates that the 5-factories can probe deep within the technically 
natural region of the 9 — ms parameter plane (see Eq. (j3J)), with sensitivity to mixing angles 
in the 10 -4 — 10~ 3 range. For light scalars with masses below the 2ir threshold we see 
that, although the branching to dimuons approaches 100%, the sensitivity is diminished 
as the S particle is very narrow and long-lived and thus able to escape the detector. We 
also observe that the sensitivity is weakened near the fo resonance, and for heavy scalars, 
as in these regions the branching to muons is small. For the axion portal, we present in 
Fig. [2] the f q — fi sensitivity for an 800 MeV pseudoscalar, indicating that the sensitivity 
to the axion couplings reaches f q j ~ 10 3 TeV. Qualitatively, we see that when f q is large, 



1 — exp 



^minFs'( a ) 

75(a) sin 9 



(13) 
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Figure 2: Sensitivity of the BaBar/Belle dataset to combined B — > Ka and B — > K*a decays in the dimuon 
channel. We have set \nA\jy/m t ~ 1. 

sensitivity is lost as the branching of B mesons to pseudoscalars is small, while for large 
// sensitivity is lost as the decays of a are primarily hadronic. Nonetheless, we note that 
the sensitivity to axion couplings obtained here appears significantly stronger than that of 
Ref. [29J. We believe that much of this numerical discrepancy can be attributed to the 
difference in experimental sensitivity to the branching fraction used in the two analyses. In 
addition, we assume at least a small hierarchy exists between the weak scale and the UV 
cutoff, whereas Ref. (29] considers the UV-complete two-Higgs-doublet model, in which - 
without this hierarchy - the top-P^ loop has an additional suppression compared to Eq. ([7]). 
Finally, it is also important to emphasize the complementarity of constraints from rare K and 
B decays. For a weakly interacting (pseudo) scalar particle with a mass below the dimuon 
threshold and a long lifetime, the K — > tc +Ip decay (e.g. K — > ttujj) is the most efficient 
probe [T7] . On the other hand, a semi-leptonic signature of S or a is more efficiently probed 
via B decays, since the CKM suppression from the top-loop is less severe. 

There are several other interesting signatures for the Higgs portal scenario in Eq. (|7|). 
Consider the decay of B° mesons to a pair of scalars. Assuming for simplicity that the A' 
trilinear vertex dominates, we obtain the following estimate for the branching to an SS pair: 



The suppression of the 2S final state relative to KS is due to the fact that the decay 
amplitude for 25* is proportional to the decay constant J'b — 200 MeV, while the KS decay 
amplitude, in the same units, is controlled by fom B ~ 2 GeV. For B d decays there is of 
course an extra CKM suppression by |V^/V^| 2 relative to (TL4"]) . Nonetheless, the overall rate 
to muons for 300 MeV scalars can reach Br^^^ ~ 10~ 8 with a moderate fine-tuning of 




(14) 
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couplings to allow for a larger A'. 

Returning to the decays mediated by XSSH^H, we note that only in the limit A > 10~ 2 
is the branching for B s — > 4fi above the 10 -8 level. Such values of A are difficult to reconcile 
with the large additive renormalization of m| by Xv 2 , which would require fine tuning at 
the level of 1 part in 10 3 for a 1 GeV scalar. Such a fine tuning can be avoided in the two- 
Higgs doublet model with a portal A/^if jifiS 12 , where Xh ± can naturally be 0(1) if tan/3 is 
maximal, tan (3 = (H2) / (Hi) ~ 50. Taking the results of the 6 — s — S 2 transition calculated 
in [H], with a charged Higgs mass = 300 GeV, we obtain the following estimate for 
the rate of the B s — > 4{i transition: 

Bv Bs ^ s ^ ~ 2 x 10- 7 x X^X'Js x Br 2 ^. (15) 

Assuming a similar sensitivity to the four-muon channel as for at CDF [39] , we conclude 
that the Tevatron experiments can probe X 2 Hl x Br|^ 2M at the 0(0.1) level. A tension in the 
parameters arises if fTTS]) is to be maximized: larger values of A^ imply larger values of ms 
where Brs_>2^ diminishes. If Bis^fi "C 1, searches for Htt + 7t~ and UK + K~ final states with 
two hadrons reconstructing the same invariant mass might be more advantageous than the 
search for fully leptonic decays of both S scalars. 

4. Rare -decays through the U(l)s sector 

In this section we will discuss 5-decays via the combined Higgs and vector portals, 

£ffiggs+Vector = -X(H ] H)(H 1 ' H') - -F^F'^, (16) 

where H' is a new scalar field charged under an additional U(l),s gauge group, while the SM is 
U(l),5-neutral. The vector portal in ffl6|) is the minimal possibility jlQ] although other options 
that involve the gauging of anomaly-free SM quantum numbers are also plausible |JT] • The 
gauging of the scalar coupled via the Higgs portal has two important consequences. First, 
as has been emphasized in many papers (see, e.g. [18]ll2ill3]), the yield of leptons in the 
final state can be enhanced, as the decay of the physical excitation h' may proceed via the 
intermediate vector states of U(l)s which in turn cascade to leptons: 

ti vv ->• nil. (17) 

The vectors decay with equal probability to different (charged) lepton species, so that the 
decay to electrons is no longer suppressed. The decay chain fTl7|) is efficient if m h i > 2m v , 
and the relative branching of V to leptons for the minimal portal is regulated by the well- 
measured process 7* — > hadrons jl3] characterized by the R{s) ratio. A second important 
consequence is that the decay chain f[T7|) is likely to be very prompt, occurring within the 
detector even for very small values of k. 

We first address B — > KV decays within the pure vector portal model. In this case, on 
account of fl2]), there is no particular enhancement. Calculation of the decay width involves 
the familiar Z and 7 penguins, with the vector particle attached via kinetic mixing. The 
result turns out to be very small, and for my ~ 1 GeV we find, 

Bt B ->kv ~ 6 x 10~V. (18) 
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This channel appears to be less sensitive to the kinetic mixing parameter k than existing 
limits from other low-energy precision experiments [T% | H5 tPPffi 

The next process we consider is B ->■ K(K*)h' -> K(K*)VV -> K(K*)llll. Utilizing the 
results (EL we obtain 



Xv'v 

m[J (l + fiK)/2) 2 ' 



^b^k{k*)uu - 0-5 x — T ) n ■ Df _ WoN9 , (19) 



having assumed that T V -^ e s = IV-^- From (fT9|) one can infer rather strong C(10 -4 ) 
sensitivity to the mixing parameter Xvv'm^ 2 . However, it is important to bear in mind that 
the naturalness limits on A are also quite strong, Xv 2 < Oim 2 ^), and therefore (|T9|) is not 
probing the natural strength of the Higgs portal. 

A particularly interesting aspect of the combined Higgs and vector portals is that the 
decay B° — > VV can proceed through an off-shell h — h' propagator. At first, it may appear 
that this process is insignificant, as both h — h' mixing and the h' — V — V vertex are 
proportional to v', naively suggesting strong suppression for a light vector. However, it 
turns out that the longitudinal vector modes in the final state cancel this v '-dependence so 
that the result remains finite in the my —> limit, 

^ * ,9,1/2 1 — 4ml/m 2 R + Y2m\rlm%> , , 

where we have taken rrih = 115 GeV. This decay leads to four leptons in the final state, and 
there is a possible enhancement of the rate for close to tub- 

Finally, the cascade decay B — > 2h' — > AV — > 4(//) leads to eight leptons in the final 
state. The rate for this process may be enhanced in the two-Higgs doublet model, and reach 
O(l0~ 7 ) x X 2 H . Therefore probes of this signature at a level better than 1 part in 10 7 at the 
Tevatron are well-motivated. 



4. Conclusions 

We have shown that rare decays of 5-mesons to semileptonic or fully leptonic final states 
can, via the 5-factory datasets, be a sensitive probe for new light states coupled through 
the Higgs and axion portals. The results of sections 2 indicate that existing data allows for 
a probe of neutral scalars coupled through the Higgs portal down to mixing angles as small 
as 10~ 3 — 10~ 4 . In addition, the axion portal coupling to the top quark can be tested at an 
impressive level of sensitivity, f q ~ 10 3 TeV. 

We have also shown that a combination of vector and Higgs portals, e.g. gauging of 
the scalar field coupled to H^H, can enhance sensitivity through the multilepton decays of 
the scalars. Among the novel signatures that we believe can be efficiently probed at both 
(super-)B factories and hadron colliders are the K(K*) +2(11), 2(11) and 4(ZZ) final states. 

4 Further model-dependent sensitivity to the kinetic mixing parameter may be obtained with cosmic- and 
gamma-ray experiments and neutrino telescopes [45] . 



9 



As far as we are aware, these final states have not been explored to date and thus represent 
a new opportunity to access light mediators. 

Finally, we should mention that while we have focused on 5-decays, and similar studies 
in the kaon sector have a long history, further sensitivity to these portal couplings may arise 
in the charm sector, via .D-decays. 

Acknowledgements 

We thank R. Kowalewski, Y. Kwon and M. Trott for helpful discussions, and especially M. 
Freytsis, Z. Ligeti and J. Thaler for emphasizing the UV-sensitivity of the b — s — a vertex in 
Eq. (7). We also thank the SLAC theory group for organizing the stimulating 'Dark Forces 
workshop' in September 2009. The work of A.R. and M.P. is supported in part by NSERC, 
Canada, and research at the Perimeter Institute is supported in part by the Government of 
Canada through NSERC and by the Province of Ontario through MEDT. 



References 

[1] B. Aubert et al [BA BAR Collaboration], Nucl. Instrum. Meth. A 479, 1 (2002) 
[arXiv:hep-ex/0105044] . 

[2] A. Abashian et al. [Belle Collaboration], Nucl. Instrum. Meth. A 479, 117 (2002). 

[3] F. Abe et al. [CDF Collaboration], Nucl. Instrum. Meth. A 271, 387 (1988). 

[4] S. Abachi et al. [DO Collaboration], Nucl. Instrum. Meth. A 338, 185 (1994). 

[5] For a recent review, see e.g. A. J. Buras, Prog. Theor. Phys. 122, 145 (2009) 
|arXiv:0904.49T7l [hep-ph] ] . 

[6] A. A. Anselm and N. G. Uraltsev, Sov. Phys. JETP 57, 1142 (1983); R. S. Willey, 
Phys. Rev. D 39, 2784 (1989); [Zh. Eksp. Teor. Fiz. 84, 1961 (1983)]; J. Prades 
and A. Pich, Phys. Lett. B 245, 117 (1990); S. N. Gninenko and N. V. Kras- 
nikov, Phys. Lett. B 427, 307 (1998); A. E. Faraggi and M. Pospelov, Phys. Lett. 
B 458, 237 (1 999); A. Pedes, H. K. Dreiner and P. Richardson, Phys. Rev. D 65, 
015001 (2001) [arXiv:hep-ph/0106199'] ; B. McElrath, Phys. Rev. D 72, 103508 (2005) 
|arXiv:hep-ph/0506151] ; P. Fayet, Phys. Rev. D 74, 054034 (2006); T. M. Aliev, 
A. S. Cornell and N. Gaur, JHEP 07 07, 072 (2007) jarXiv: 0705 .45421 [hep-ph]]. 

[7] S. Chang, R. Dermisek, J. F. Gunion and N. Weiner, Ann. Rev. Nucl. Part. Sci. 58, 75 
(2008) jarXiv:0801.4554l [hep-ph]]. 



[8] B. Patt and F. Wilczek, arXiv:hep-ph/0605188 



[9] R. Foot, H. Lew and R. R. Volkas, Phys. Lett. B 272, 67 (1991). R. Foot and X. G. He, 
Phys. Lett. B 267, 509 ( 1991); D. G. Cerdeno , A. Dedes and T. E. J. Underwood, 
JHEP 0609, 067 (2006) |arXiv:hep-ph/0607157]; J. R. Espinosa and M. Quiros, Phys. 



10 



Rev. D 76, 076004 (2007) |arXiv:hep-ph/0701145| ; J. March-Russell, S. M. West, 
D. Cumberbatch and D. Hooper, JHEP 0807, 058 (2008) jarXiv:0801.3440l [hep-ph]]; 
M. Ahlers, J. Jaeckel, J. Redondo and A. Ringwald, Phys. Rev. D 78, 075005 (2008) 
|arXiv:0807.4T43l [hep-ph]]; J. L. Feng, H. Tu and H. B. Yu, JCAP 0810, 043 (2008) 
|arXiv:0808.2"3T8l [hep-ph]]; K. Kohri, J. McDonald and N. Sahu, larXiv:0905.13T2l 
[hep-ph]; J. L. Feng, M. Kaplinghat, H. Tu and H. B. Yu, JCAP 0907, 004 (2009) 
jarXiv:0905.3039l [hep-ph]]. 

[10] O. Adriani et al. [PAMELA Collaboration], Nature 458, 607 (2009) [ arXiv:0810.4995l 
[astro-ph]]. 

[11] N. Arkani-Hamed, D. P. Finkbeiner, T. R. Slatyer and N. Weiner, Phys. Rev. D 79, 
015014 (2009); M. Pospelov and A. Ritz, Phys. Lett. B 671, 391 (2009) [arXiv:0810.1502l 
[hep-ph]]; Y. Nomura and J. Thaler, Phys. Rev. D 79, 075008 (2009) [arXiv:0810.5397l 
[hep-ph]]. 

[12] R. S. Willey and H. L. Yu, Phys. Rev. D 26, 3287 (1982). 

[13] B. Grinstein, L. J. Hall and L. Randall, Phys. Lett. B 211, 363 (1988). 

[14] C. Bir d, P. Jackson, R. V. Ko walewski and M. Pospelov, Phys. Rev. Lett. 93, 201803 
(2004) [arXiv:hep-ph/0401195|; C. Bird, R. V . Kowalewski and M. Pospelov, Mod. Phys. 
Lett. A 21, 457 (2006) [arXiv:hep-ph/0601090j ; G. K. Yeghiyan, larXiv:0909.4"9l9l [hep- 
ph]; A. Badin, G. K. Yeghiyan and A. A. Petrov. larXiv:0909.5219 l [hep-ph]: C. S. Kim, 
S. C. Park, K. Wang and G. Zhu, larXiv:0910.429T1 [hep-ph]. 

[15] D. O'Connell, M. J. R amsey-Musolf and M. B. Wise, Phys. Rev. D 75, 037701 (2007) 
[arXiv:hep-ph/0611014| . 

[16] G. Hiller, Phys. Rev. D 70, 034018 (2004) [ arXiv:hep-ph/0404220] ; Z. Heng, R. J. Oakes, 
W. Wang, Z. Xiong and J. M. Yang, Phys. Rev. D 77, 095012 (2008) [ arXiv: 0801 .11691 
[hep-ph]]. 

[17] M. Pospelov, A. Ritz and M. B. Voloshin, Phys. Lett. B 662, 53 (2008) [arXiv:0711.4"866l 
[hep-ph]]. 

[18] M. Pospelov, Phys. Rev. D 80, 095002 (2009) [arXiv: 081 1.10301 [hep-ph]]. 

[19] M. J. Strassler and K. M. Zurek, Phys. Lett. B 651, 374 (2007) |arXiv:hep-ph/0604261| . 

[20] N. Arkani-Hamed and N. Weiner, JHEP 0812, 104 (2008) [ arXiv:0810.07T4l [hep-ph]]. 

[21] V. Silveira a nd A. Zee, Phys. Lett. B 161, 136 (1985); J. McDonald, Phys. Rev. D 50, 
3637 (1994) [arXiv:hep-ph/0702143] . 

[22] C. P. Burgess, M. Po spelov and T. ter Veldhuis, Nucl. Phys. B 619, 709 (2001) 
[arXiv:hep-ph/0011335l ; 



11 



[23] G. Anderson and L. Hall, Ph ys. Rev . D 45, 2685 (1992); J. Espinosa and M. Quiros, 
Phys. Lett. B 305, 98 (1993) |arXiv:hep-ph/9301285] ; S. Profumo, M. Ramsey-Musolf 
and G. Shaughnessy, JHEP 0708, 010 (2007) jarXiv: 0705 .24251 [hep-ph]]; T. Clark, 
B. Liu, S. Love and T. ter Veldhuis, Phys. Rev. D 80, 075019 (2009) [ arXiv:0906.5595l 
[hep-ph]]; R. Lerner and J. McDonald. arXiv:0909.0520l [hep-ph]. 

[24] N. V. Krasniko v, Phys. Lett. B 2 91, 89 (1992); N. V. Krasnikov, Mod. Phys. Lett. A 
13, 893 (1998) |arXiv:hep-ph/9709467|. 



O. J. P. Eboli and D. Zeppenfeld, Phys. Lett. B 495, 147 (2000) |arXiv:hep-p h/0009158|. 



R. Schabinger and J. D. Wells, Phys. Rev. D 72, 093007 (2005) |arXiv:hep-ph/0509209|. 
M. Bowen, Y. Cui and J. D. Wells, JHEP 0703, 036 (2007) [arXiv:hep-ph/0701035X 



M. J. Strassler and K. M. Zurek, Phys. Lett. B 661, 263 (2008) [arXiv:hep-ph/0605193] . 

V. Barger, P. Langacker, M. McCaskey, M. J. Ramsey-Musolf and G. Shaughnessy, 
Phys. Rev. D 77, 035005 (2008) [ arXiv:0706.4"3lT1 [hep-ph]]; V. Barger, P. Langacker, 
M. McCaskey, M. Ramsey-Musolf and G. Shaughnessy, Phys. Rev. D 79, 015018 (2009) 
jarXiv:0811.0393l [hep-ph]]. 

M. Freytsis, Z. Ligeti and J. Thaler, Phys. Rev. D 81, 034001 (2010) |arXiv:0911.5"355l 
[hep-ph]]. 

A. Ali, P. Ball, L. T Handoko and G. Hiller, Phys. Rev. D 61, 074024 (2000) 
|arXiv:hep-ph/991022l"] . 



P. Ball and R. Zwicky, Phys. Rev. D 71, 014015 (2005) | arXiv:hep-ph/04 06232|; P. Ball 
and R. Zwicky, Phys. Rev. D 71, 014029 (2005) larXiv:hep-ph/0412079] . " 

C. Amsler et al. (Particle Data Group), Phys. Lett. B 667, 1 (2008). 

B. Aubert et al. [B ABAR Collaboration], Phys. Rev. D 73, 092001 (2006) 
|arXiv:hep-ex/0604007] . 

J. T. Wei et al. [BELLE Collaboration], Phys. Rev. Lett. 103, 171801 (2009) 
jarXiv:0904.0770l [hep-ex] ] . 

H. J. Hyun et al. [Belle Collaboration], Phys. Rev. Lett. 105, 091801 (2010) 
[arXiv:1005. 1450l [hep-ex]]; Y. Kwon, talk at the Dark Forces Workshop, SLAC, Sept 
24-26, 2009, http ://indico . cern . ch/ conf erenceDisplay. py ? conf!d= 6 7760 1 

M. B. Voloshin, Sov. J. Nucl. Phys. 44, 478 (1986) [Yad. Fiz. 44, 738 (1986)]. 

J. F. Gunion, H. E. Haber, G. Kane and S. Dawson, The Higgs hunter's guide (Addison- 
Wesley, 1990). 

[38] T. N. Truong and R. S. Willey, Phys. Rev. D 40, 3635 (1989). 



12 



[39] T. Aaltonen et al [CDF Collaboration], Phys. Rev. Lett. 100, 101802 (2008) 
|arXiv: 0712 .17081 [hep-ex]]. 



[40] B. Holdom, Phys. Lett. B 166, 196 (1986). 

[41] L. B. Okun, Sov. Phys. JETP 56, 502 (1982) [Zh. Eksp. Teor. Fiz. 83, 892 (1982)]; 
R. Foot, X. G. He, H. Lew and R. R. Volkas, Phys. Rev. D 50, 4571 (1994) 



arXiv:hep-ph/9401250|; P. J. Fox and E. Poppitz, Phys. Rev. D 79, 083528 (2009) 



arX iv:0811 .0399 [hep-ph]]. 



[42] S. Gopalakrishna, S. Jung and J. D. Wells, Phys. Rev. D 78, 055002 (2008) 
|arXiv:0801.3456l [hep-ph]]. 

[43] B. Batell, M. Pospelov and A. Ritz, Phys. Rev. D 79, 115008 (2009) [ arXiv: 0903 .03631 
[hep-ph]]; R. Essig, P. Schuster and N. Toro, Phys. Rev. D 80, 015003 (2009) 
[ar Xiv:0903.394H [hep-ph]]; M. Reece and L. T. Wang, JHEP 0907, 051 (2009) 
|arXiv:0904. 17431 [hep-ph]]. 

[44] J. D. Bjorken, R. Essig, P. Schuster and N. Toro, Phys. Rev. D 80, 075018 (2009) 
|arXiv:0906.0580l [hep-ph]]; B. Batell, M. Pospelov and A. Ritz, larXiv:0906.5"6T4l [hep- 
ph]. 

[45] B. Batell, M. Pospelov, A. Ritz and Y. Shang, larXiv:0910. 15671 [hep-ph]; P. Schuster, 
N. Toro and I. Yavin, arXiv:0910.1602 [hep-ph]; P. Schuster, N. Toro, N. Weiner and 
I. Yavin. larXiv :0910. 1839 [hep-ph]; P. Meade, S. Nussinov, M. Papucci and T. Volansky, 
larXiv:0910.4"T60l [hep-ph] . 



13 



